Measurement of BB Angular Correlations 
at ^/s = 7 TeV with the CMS Experiment 
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Measurements of the angular correlations between beauty and anti-beauty hadrons produced 
in LHC pp collisions at ^/s = 7 TeV are presented. These results probe for the first time 
the small angular separation region and show sensitivity to coUinear particle emission. The 
results are compared with predictions based on perturbative QCD calculations at leading and 



X 

• next-to-leading order. 

1 Introduction 

l/^ I Studies of production properties of beauty quarks (b) at the CERN LHC collider are of 

I twofold interest. Firstly, the b production process provides an excellent opportunity to 

study details of perturbative Quantum Chromodynamics (pQCD). Over the years, the var- 
ious tensions between the predictions and the measurements, that existed in data at lower 
energies such as the HERA or the Tevatron collider, have been reduced, however not com- 
pletely resolved. Studies at the LHC collider with higher centre-of-mass energies comple- 
ment the previous data, but also expand the reach and provide tests at precisions below 
^ ■ the present theoretical uncertainties. 

H ■ 

I Secondly, b quark production constitutes one of the major backgrounds in many of the 

searches for new physics. Any production channel of exotic states, that produces top 
quarks or W-bosons, will inherently have a large b production rate. It is of importance 
not only to understand the absolute production rates, but also to be able to describe the de- 
tails of the b production dynamics. Thus, a solid understanding of the topology of the final 
states will be crucial to constitute efficient criteria to distinguish possible signal signatures 
from b-induced background configurations. 

Within the leading-order (LO) QCD picture, the production of bb in pp collisions at LHC 
can be attributed to three parton level production subprocesses, commonly denoted by 
flavour creation (FCR), flavour excitation (FEX) and gluon splitting (GSP) (see Fig. 1). At 
higher orders, the distinction becomes scale dependent, and is thus less well defined. Due 
to the different dynamics of these components, the final state topologies differ substantially 
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Figure 1 : Examples of schematic Feynman diagrams for the three subprocesses: 
flavour creation (FCR, top), flavour excitation (FCR, middle) and gluon splitting 
(GSP, bottom). 



from each other. FCR pertains to the 2 — )■ 2 processes gluon-fusion and qq annihilation, 
where the b and b are emitted in a back-to-back configuration. FEX refers to the 2 — > 
3 process, where one b quark of a bb pair from the proton sea participates in the hard 
scattering, thereby producing an asymmetry in the momentum and angular distribution 
of the final state. The GSP contribution on the other hand describes gluons from either 
initial or final state, that split into a bb pair, which in turn are emitted preferentially at 
small opening angles and low px- Furthermore, the relative production rates themselves 
vary also as a function of the energy scale. It is expected, that at higher energies the gluon 
splitting contributions dominate, ie. processes with a coUinear branching of gluons into 
BB pairs will become the major source of b quark production. 



2 Measurements of 66 Angular Correlations 

CMS has performed the first measurement [1] of the angular correlations between beauty 
and anti-beauty hadrons (BB) produced in pp collisions at i/s = 7 TeV, thereby probing for 
the first time the region of small angular separation. The analysis is based on a data sample 
corresponding to an integrated luminosity of 3.1 ± 0.3pb^^. A detailed description of the 
CMS detector can be found in Ref. [2]. 

The measurements are done differentially as a function of the opening angle for different 
event scales, which are characterised by the leading jet transverse momentum in the event 
(independently of b hadrons). The leading jet of the event is used to trigger. The trigger 
thresholds are chosen such as to reach an efficiency over 99% for all three energy scale bins, 
which correspond to a leading jet pj in excess of 56, 84 and 120 GeV, respectively, when 
using corrected jet energies. 

The cross sections are determined by applying efficiency corrections and normalising to 
the total integrated luminosity. The angular correlations between the two B hadrons are 
measured in terms of the difference in azimuthal angles (A<^) in radians and the combined 
separation variable AR = a/ (Af/)-^ + (A^)^, where A;/ is the pseudorapidity. The analy- 
sis results are quoted for the visible kinematic range defined by the phase space at the 
B hadron level by the requirements |j/(B)| < 2.0 and Pt(B) > 15 GeV for both of the B 
hadrons. The leading jet used to define the minimum energy scale is required to be within 
a pseudorapidity of |?/(jet) | < 3.0. 
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CMS \ls = 7 TeV, L = 3.1 pb^ 
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Figure 2: Differential BB production cross sections as a function of (left) and lS.(p (rad) 
(right) for the three leading jet regions. For clarity, the pj > 56 and 84 GeV bins are 
offset by a factor 4 and 2, respectively. The error bars of the data show the statistical 
(inner) and the total (outer bars) uncertainties. A common uncertainty of 47% due to 
the absolute normalisation on the data points is not included. The PYTHIA prediction is 
normalised to the region AR > 2.4 or /^(p > 2.4. The widths of the shaded bands indicate 
the statistical uncertainties of the predictions. 



In order to measure the angular correlations also in the collinear regime, the reconstruction 
of the B hadrons is done independently of jet algorithms. The method uses the B hadron 
decays and is based on an iterative inclusive secondary vertex finder that exploits the ex- 
cellent CMS tracking information [3]. This allowed to approximate the flight direction of 
the original B hadron by the vector between the primary (PV) and the secondary vertices 
(SV). A resolution of 0.02 rad in AR could be achieved that way. The average overall event 
reconstruction efficiencies (for both B hadrons) are found to be of order 10% at an average 
purity of 84%. Detailed studies were performed to ensure high accuracy in the B-hadron 
kinematics description. In addition, the angular dependence of the efficiency description 
was verified by a special event mixing technique, both in data and the simulation. 

The measured cross sections are presented in Fig.2. Overlaid are the predictions by the 
PYTHIA calculations, which are normalised to the AR > 2.4 or A^ > 2.4 (rad) regions, 
where the calculations are expected to be more reliable. Note, that an overall common 
uncertainty of 47% due to the absolute normalisation is not shown in the figures. 

We find that the cross sections at small AK or Acp are substantial and even exceed the 
values observed at large angular separation values. Hence, the configurations where the 
two B hadrons are emitted in opposite directions are much less likely than the collinear 
configuration. 

The measurements are compared to various predictions, based on LO and next-to-leading 
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Figure 3: Ratio of the differential BB production cross sections, as a function of AR (left) 
and (rad) (right), for data, MadGraph, MC@NL0 and CASCADE, with respect to the 
PYTHIA predictions, shown also for the three leading jet pj bins. The simulation is nor- 
malised to the region AR > 2.4 and A(^ > 2.4 (rad) (FCR region), as indicated by the 
shaded normalisation region. The widths of the theory bands indicate the statistical un- 
certainties of the simulation. 



(NLO) pQCD calculations. Figure 3 illustrates the shape sensitivity by showing the ratio of 
the different AR distributions to the PYTHIA Monte Carlo predictions. It is found, that the 
overall tendency in shape is in general reasonably described by the predictions, however 
the normalizations and the details in shape, in particular at small opening angles are not 
described well by any of the calculations. Apart from MadGraph program, all predic- 
tions underestimate the amount of gluon splitting contributions in the collinear region. 

Perturbative QCD predicts a back-to-back configuration for the production of the BB pair 
(i.e. large values of AR and/or A^) for the LO processes. In contrast, the region of phase 
space with small opening angles between the B and B hadrons provides strong sensitivity 
to collinear emission processes, such as the ones present in higher-order processes. Gluon 
radiation which splits into bb pairs is anticipated to have a smaller angular separation 
between the b quarks. 

The measurements show that the BB production cross section ratio p/sR = crAR<o.8 / '^ar>2.4 
increases as a function of the leading jet pj in the event (see Fig. 4). Larger px values lead 
to more gluon radiation and, hence, are expected to produce more gluon splitting into 
BB pairs. This general trend is described reasonably by the theoretical calculations. 
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Figure 4: Left: ratio between the BB production cross sections in AR < 0.8 and AR > 2.4, 
P^R = c^AR<0.8 / o^^R>lA, as a function of the leading jet pj. Right: asymmetry between 
the two regions, (crAR<o.8 - c^ar>2.4) / (c^ar<0.8 + '^ar>ia)- The symbols denote the data 
averaged over the bins and are plotted at the mean leading jet px of the bins. For the data 
points, the error bars show the statistical (inner bars) and the total (outer bars) errors. 
Also shown are the predictions from the PYTHIA and MadGraph simulations, where 
the widths of the bands indicate the imcertainties arising from the limited number of 
simulated events. 

3 Conclusions 

The first measurements of inclusive beauty production have been performed at the LHC by 
the CMS experiment over a large range from very low transverse momenta up to 300 GeV 
in the central rapidity region. Comparisons with theoretical predictions, based on pQCD 
calculations have confirmed the large production cross section. The calculations in general 
describe the overall features of beauty production fairly well. However, the predictions 
do not yet adequately describe the differential measurements, neither in the B transverse 
momentum, nor the rapidity nor the BB opening angle distributions. 
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